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Abstract. Wc present a simple approach to study the one-dimensional pressureless Euler 
system via adhesion dynamics in the Wasserstein space CP2(M) of probabihty measures with 
finite quadratic moments. 

Starting from a discrete system of a finite number of "sticky" particles, we obtain new 
explicit estimates of the solution in terms of the initial mass and momentum and we are able to 
construct an evolution semigroup in a measure-theoretic phase space, allowing mass distributions 
in CP2(IR) and corresponding L^-velocity fields. Wc investigate various interesting properties of 
this semigroup, in particular its link with the gradient flow of the (opposite) squared Wasserstein 
distance. 

Our arguments rely on an equivalent formulation of the evolution as a gradient flow in the 
convex cone of nondecreasing functions in the Hilbert space L'^(0, 1), which corresponds to 
the Lagrangian system of coordinates given by the canonical monotone rearrangement of the 
measures. 



1. Introduction 

In the recent years considerable attention has been devoted to the 1-dimensional pressureless 
Euler system 

dtp + dx{pv) = 0, 
dt{pv)+d,{pv^)^Q, 



(1-1) So/ N , a / ,2x n inMx(0,+oo); = po, '^\t^Q='"o, 



in connection with the Zeldovich model [30] for the evolution of a "sticky particle system" (SPS, 
in the following) via adhesion dynamics. This model describes the behaviour of a finite collection 
of particles, freely moving in absence of forces and sticking under collision; they can be mathe- 
matically represented by a time-dependent discrete measure p^ :— X^ILi ''^i^xi{t) concentrated in 
a finite set of 

N particles Pi{t) :~ (mi, Xi{t),Vi{t)), i = 1, . . . , N, with positive mass jti,;, ordered 
positions xi{t) < X2{t) < ... < XN-iit) < XAr(t), and velocities Vi{t). 
Denoting by Ji{t) {j : Xjit) = Xi(t)} the collection of (the indexes of) the particles Pj{t) 
coinciding with Pi{t) at the time t, the adhesion dynamic imposes that the sets Ji{t) are non- 
decreasing in time, so that Wj(t+) = ^'i(^+) for every j S Ji{t). We can thus order in a finite 
and monotone sequence < ti < t2 < ■ ■ ■ the collection of times when the cardinality of some 
Ji{t) has a discontinuity (corresponding to some collision). In each open interval [tfc,tfc+i) the 
(right-continuous) velocities Vi{t) = Xi{t) are thus supposed to be constant, and at each collision 
time tk the conservation of mass and momentum yields the update equation for the velocities 

(1.2) v.it,+)= ^%^'^^''''^^^''-\ . = l,...iV. 

It is not difficult to check that the measures and (pv)^ X^iLi '^a;i(t) solve (1.1). 

Starting from the discrete SPS, existence of measure valued solutions to (1.1) with general initial 
data and satisfying suitable entropy conditions [5, Bouchut] has been proved by Grenier [14] 



Key words and phrases. Pressureless Euler equation. Sticky particles, Wasserstein distance. Monotone rearrange- 
ment. Gradient fiows. 

G.S. has been partially supported by MIUR-PRIN'06 grant for the project "Variational methods in optimal 
mass transportation and in geometric measure theory" . 



2 



LUCA NATILE AND GIUSEPPE SAVARE 



and E, Rykov & Sinai [28] (but sec also the contribution of Martin &; Piasecki [16]) as limits 
(in the sense of weak convergence of measures) of the discrete particle evolutions as TV t +00. 
Here we also quote the different approaches of Bouchut & James [6], of POUPAUD & Rascle 
[19], and of Sever [24] in the multidimensional case; viscous regularizations of (1.1) have been 
studied by Sobolevskii [26] and Boudin [7], and a different model, starting from particles of 
finite size, has been considered by Wolansky [29] . 

The convergence result has further been extended and refined by Brenier & Grenier [9], 
Huang & Wang [15], and Nguyen & Tudorascu [18] (by a different probabilistic approach 
MOUTSINGA [17] has recently been able to consider initial velocities with nonpositive jumps at 
each points of the support of po): the basic assumption is that the discrete initial velocity Vi is the 
value in Xi of a given continuous function v with at most linear growth, and (the total mass being 
normalized to 1) the sequence converges to po w.r.t. the L^-Wasserstein distance in the space 
J'2(K) of probability measures with finite quadratic moment. This includes the case (considered in 
[9] ) of a sequence with uniformly bounded support and weakly converging to pa in the duality 
with continuous real functions. 

All these results depend on a remarkable characterization of the solution p found by Brenier & 
Grenier [9]: by introducing the cumulative distribution function Mp associated to a probability 
measure p £ 3'(R) 

(1.3) Mp{x) := p{{- 00, x\) e M, sothaXp^d^Mp in ^'(M), 

they prove that the function M{t, ■) AIp^(-) is the unique entropy solution of the scalar conser- 
vation law 

(1.4) dtM + dtA{M) = in M x (0, +00), 

where ^ : [0, 1] ^ M is a continuous flux function depending only on po and vq (see Theorem 6.1 
for a precise statement). 

It can also be shown [18] that this solution satisfies the Oleinik entropy condition 

(1.5) Vt{x2) - vt{xi) < j{x2 - xi) for pt-a.e. Xi,X2 eM., xi < X2- 

In the present paper we discuss various refinement of Brenier-Grenier result by a different 
approach. Our starting point (Theorem 2.2) is an explicit Lipschitz estimate (in the iP-Wasserstein 
distance Wp for every p > 1, see (2.1)) of the dependence of pt with respect to the initial data 
Po, {pv)q- for p = 2 it shows that {p^)^^-^ is a Cauchy sequence in T2(M) and in particular yields 
the convergence results of [28, 9, 18] allowing general initial measures in J'2(R) and (possibly 
discontinuous) velocity field vq S L'^{po). We also show that a suitable X^-likc integral distance 
between the momentum pv of two solutions can be controlled in terms of the initial data and prove 
further precise representation properties of the solution and its velocity field (Theorem 2.3). 

This leads to the construction of a semigroup S't associated to the evolution of SPS, which ex- 
hibits interesting links with another semigroup (recently studied by Ambrosio, Gigli & Savare 
[1] ) , obtained as the gradient flow in J'2 (R) of the (opposite) squared Wasscrstein distance from a 
fixed reference measure. 

This link (which at a first sight may look unexpected) can be better understood in the simpler 
case when the initial velocity field v satisfies a one-sided monotonicity condition (see section 5.4.2 
of Villani's book [27] for more details): still considering the simpler discrete case, if 

(1.6) -<5-i := min ^(^'^ ~ ^(^^^ < 0, v{x^) := Vi, 

for t G [0,S) the map Xq(x) := x + tv{x) is nondecreasing on the support of po (the finite set 
{xi : i = 1, . . . N}), so that the first collision occurs at t := S and in the interval [0, S) one has the 
freely moving measures 

N N 

(1-7) Pt := {xD^pq ^^rUiS^.+tv,, {pv)t = miViSx+tv, , t e[Q,5), 

i=l i=l 
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solving the prcssurclcss Euler system (1.1). On the other hand, the curve t ^ pt, t E [0,6], is 
a constant speed minimal geodesic in ^'2(1^) connecting po with rj := ps', as in any Riemannian 
manifold, it coincides (up to a suitable rescaling, [1, Theorem 11.2.10]) with the gradient flow in 
J'2(]K) of the functional (f)''"{p) := — ^^^(p, p°). After the collision at time t = 5 the trajectory 
of the gradient flow does not coincide with the free motion (1.7) anymore, since its velocity has a 
jump which can be described exactly by (1.2) [1, Theorem 10.4.12]. At a later time, the velocity 
field induced by the (rescaled) Wasserstein gradient flow can be characterized by the formula 

(1.8) v.{t+) = t-^ Ut) , ^^l^...N, 

and it is an interesting property, stated in Theorems 2.4 and 2.5, that the two different laws (1.2) 
and (6.3) give rise to the same evolution, even for arbitrary initial data. 

In order to obtain these results, we adopt the point of view of 1-dimensional optimal trans- 
portation and we represent each probability measures p S 3'2 (I^) by their monotone rearrangement 
Xp, which is the pseudo-inverse of the distribution function Mp of (1.3) (a similar approach, in a 
probabilistic framework, has been also used by [17]; see also [13] for other applications) 

(1.9) Xp{w) ■.= mf{x : Mp{x) > w} =M{x : p{{-oo,x]) > w} we (0,1). 

The map p 1-^ Xp is an isometry between 5'2(K) (endowed with the i^-Wasserstein distance) 
and the convex cone % of nondecreasing functions in the Hilbert space L^{0, 1). Through this 
isometry, any gradient flow with respect to W2 in 5'2(If^) can be rephrased as a gradient flow in 
3C with respect to the L^{Q, l)-distance, and one can use the powerful tools the classical theory of 
variational evolution inequalities in Hilbert spaces (we refer to the book by Brezis [10]). It turns 
out (see Theorem 2.6) that in this Lagrangian formulation the solution Xp^ admits three simple 
characterizations, in terms of the L^{0, l)-projection Pjc onto X 

(1.10) Xp^ = PjciXp, + tVo), Vo^voo Xp„ 
and of the differential inclusions 

I-x being the indicator function associated to % (see next (2.34)). (1.10) and (1.11) encode all 
the qualitative information on the measure- valued solution pt, and their proof in the case of the 
discrete SPS constitutes the core of our argument. It relies on an elementary but careful description 
of the L^(0, l)-projection operator Pjc and on the subdifferential of I-x, which has been carried 
out in Section 3. Once pt has been determined, its velocity vt £ ip,(R) can be recovered from 

the right derivative V{t) := ^Xp^ € i^(0, 1): in fact, as a byproduct of the second differential 
inclusion of (1-11), V{t) is a function of X{t) and therefore one obtains 

(1.12) Vit)^vtoXp^. 

The projection formula (1.10) (which has been introduced by Shnirelman [25, 2] in a slightly 
different form, see Remark 2.9) lies more or less explicitly at the core of the formulations by [28] and 
[9]. As it has been nicely explained by Andrievsky, Gurbatov & Sobolevskii [2] elaborating 
the contribution of [25], (1.10) is equivalent to the Generalized Variational Principle of [28], which 
can be expressed through the convex envelope of the primitive function of the map + tVo: as 
stated in full generality by Theorem 3.1, this convexiflcation characterizes the L^-projection on OC. 
On the other hand, a convexiflcation is also involved in the second Hopf formula for the solutions 
of the Hamilton- Jacobi equation associated to (1.4), as it has been already observed by [9, §4]: we 
will detail this point in Theorem 6.1. 

The link between the formulation based on the scalar conservation law (1.4) and the Hilbertian 
theory of gradient flows like (1.11) is not at all surprising, after the illuminating paper by Bre- 
nier [8]. Wasserstein contraction properties of solutions of one-dimensional scalar conservation 
laws have also been recently obtained by Bolley, Brenier & Loeper [4] (see also the further 
contribution by Carrillo, Di Francesco & Lattanzio [11]). So it would be possible in prin- 
ciple to approach SPS starting from (1.4) and trying to apply the techniques developed there. 



(1-11) :r.^p. + dWp,) 3 Vo, t-Xp, + diociXp,) 3 Xp, - Xp„, 
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Notice however that two solutions originating from different initial distributions of position and 
velocity give rise to two scalar conservation laws differing not only by the initial data but also 
by the flux functions, so that their comparison does not look immediate. Moreover, the present 
self-contained approach is very simple, since it relies on elementary tools of convex analysis and 
direct computations on the discrete case; the simultaneous characterization of the evolution by 
(1.10) and (1.11) provides a more refined description of the solution and, as a byproduct, a new 
direct proof of Brenier & Grenier theorem. 

Plan of the paper. In the next section we recall some basic definition and notation and we 
state our main results. Section 3 collects the main properties related to the convex cone OC in 
L^(0, 1) (projection, polar cone, subdifferential of the indicator function): they provide simple 
but crucial tools for the analysis of the discrete SPS presented in Section 4, which contains all 
the basic calculations. Section 5 deals with existence, stability, and uniqueness of the solution in 
the Lagrangian formulation. The final steps of the proofs (mainly concerning the various limit 
processes) will be detailed in the last Section 6, where we also show a new derivation of Brenier 
&; Grenier Theorem [9] from the Lagrangian representation of the SPS. 

2. Main results 

Couplings, Wasserstein distance, and monotone rearrangementa. For p G [l,+cx3) let us 

denote by J'p(M) the space of Borel probability measures p with finite p-moment \x\p dp{x) < 
+00. The LP Kantorovich-Rubinstein- Wasserstein distance Wp{p^,p^) between two measures 
p^,p^ e 5'p(K) can be defined in terms of couplings, i.e. probability measures p G ^(R x R) 
such that TT^p = (0*, I = 1, 2, by the formula 

(2.1) WPip\p'):=mm{ f |x - dp(a;, y) : p G T(R x R), Tr^p = p'}. 

Here n'^{xi,X2) = Xi is the usual projection on the i-th coordinate and for a general Borel map 
T : R™ R" and a Borel measure p G CP(R'") the push-forward v = T^p is the measure defined 
by viA) = p{T^^{A)) for every Borel set A C R". We will repeatedly use the change-of-variable 
formula 

(2.2) J C{y) d{T#p){y) = J C{T{x)) dp{x) for every Borel map C : R" ^ [0, +oo]. 

More generally, given a convex, even, and lower semicontinuous function ?/; : R ^ [0, +c»], we can 
consider the cost c^{x,y) := iplx — y), a;,y G R, and the associated optimal mass transportation 
problem 

(2.3) e,^(p\p2) :^inf| / xb{x - y) dp{x,y) : p e Ti^ x R), n^^p ^ p^} . 

In the present 1-dimensional case, there exists a unique optimal coupling p = To{p^,p^) realizing 
the minimum of (2.1) and of (2.3) (at least when the cost is finite): it can be explicitly characterized 
by inverting the distribution functions of p^,p^. More precisely, for every p G J'(R) we consider 
its monotone rearrangement Xp (1.9), a right-continuous and nondecreasing function satisfying 

(2.4) {Xp)#X = p, A:=^i|^^^^^, £c{x)dp{x)= l\{Xp{w))dw 

for every nonnegative Borel map C : R ^ [0, +0x3]. In particular, p G 5'p(R) iff Xp G LP{0,1). 
Moreover, thanks to the Hoeffding-Frechet theorem [20, Sec. 3.1], the joint map Xpi_p2{w) := 
{X pi (w) , X p2 (w)) , w G (0,1), characterizes the optimal coupling p G To{p^,p^) by the formula 

(2.5) p={Xpi^p.)^X, 
so that [12, 20, 27] 

(2.6) WP{p\p^)^ [ \Xpiiw)-Xp2{w)\''dw, eip\p^)= [ i^{Xpi{w)-Xp2{w))dw, 
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and the map p E CP(]R) i — > Xp is an isonietry between J'2(I^) and the convex subset X of L^{0, 1) 
of (essentiahy) nondeereasing functions (which can be identified with their right-continuous rep- 
resentatives) . 

An explicit estimate through Wasserstein distance. We introduce the set 

(2.7) Vp(R) := (p,pi;) G J'p(M) X M(K) : w e LP(R)}, pe[l,+oo), 

M(M) being the set of aU signed Borel measures with finite total variation, the semi-distances 
(here = (p*, p*w*)) 

(2.8) m{n\p?):^[ \v\x)-v\yWdp{x,y) p^T,{p\p^) 



(2.9) =/ \v\Xp.{w))~v\Xp.{wWAw, 

Jo 

and the distances 

(2.10) DPip\p') :=W^;(p\ p^) + U^{fi\p'). 
We also set 

(2.11) [p];:= [ {\x\P + vP{x))Ap{x)=Dl{t,,{5,^Q)). 

Proposition 2.1. Dp is a distance in Vp(]R) and (VplM), Dp) is metric (but not complete) space 
whose topology is stronger than the one induced by the weak convergence of measures. The collec- 
tion of discrete measures 

N N N 

(2.12) V(R) := |/i= (^m,(5^,,^m,w,(5^J : > 0, ^m, = 1, Xi,iH€Rj 

i—l i—1 i—1 

is a dense subset o/Vp(R). A sequence pn = (Pn, PnVn), n €N, converges to p = (p, pv) in Vp(M.), 
p > 1, if and only if (see [1, Def. 5.4.3] j 



(2.13) H^p(p„,p)^0, pnVn^pv weakly in M{R), / |w„|^dp„^ / |w|^dp. 

Jr Jr 

Let us denote by S^t ■ V(M) V(R) the map associating to any discrete initial datum (po, PqVq) 
the solution [pt,ptVt) of the (discrete) sticky-particle system. S^t is a semigroup in 



Theorem 2.2 (Stability with respect to the initial data). Let pi = {pl^pfvf) = 5^t[Po\, ^ = 1,2, 
be the solutions of the (discrete) sticky-particle system with initial data p^ 6 V(R). Then for every 
convex cost (2.3) and every p >1 



(2.14a) e^(pi,p?)</ i:{x + tv\x)-{y + tv''{y))&p{x,y), p = Toip\p^), 

jRxR 

(2.14b) Wpiplpj) < Wp{pl,pl)+tUp{pl,pl), 

(2.14c) ^ Ul{pl,pl)dr < C{l + t)(^[p% + [p%)(w2{pI,pI) + U2{pI,pI)), 

for a suitable "universal" constant C independent of t and the data. 

We say that a map : Vp(R) — » Vp(R) is strongly-weakly (s-w) continuous if for every p^ , p E 
Vp(R) with ^[p"] = (p",p"w"), .y[p] = ip,pd) G Vp(R), 

(2.15) lim Dp{pn, p) = =^ lim Wp(p",p) = 0, PnVn ^ pv weakly in 

nf-t-oc nf-t-oo 

Theorem 2.3 (The evolution semigroup in Vp 
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(a) The semigroup S^t can he uniquely extended by density to a right- continuous semigroup (still 
denoted S^t) of strongly-weakly continuous transformations in Vp(M), p >2, thus satisfying 

(2.16) S^s+M^^A'^tM Vs,t>0, \imDp{yt[lA,l^) = VMeVp(M). 

complies with the same estimates (2.14a, 6, c) of Theorem 2.2. 
(h) (pt, pfVt) = S^titA^ M G "^20^), "is a distributional solution of (1.1) satisfying Oleinik entropy 
condition (1.5). 

(c) If 4' : W R is a convex function such that ipivo) € L^^{M.), and (pt^ptVf) = S^Ai^q], then 

(2.17) the map t t-^ ip{vt) dpt{x) is nonincreasing in [0, +oo), 



and its (at most countable) jump set 7 = T(/i) is independent ofijj. 

(d) If pi € Vp(IR) and pt = {pt,PtVt) = -^t[A*], t G [0, +oo), the curve t ^ pt is Lipschitz in CPp(R) 
with respect to Wp, and the curve t ^ ptVt is continuous with respect to the weak topology in 
M(M), right- continuous in [0, +c») with respect to the (semi-) distance Up, and left- continuous 
at each t £ (0, +oo) \ T where 1 is the at most countable jump set of (2.17). 

(e) Let p" = {p",p"v") = S^tip"] and pt = {pt,ptVt) = 'S^tllA'j */ converges to p in Vp(R) as 
n t +00, then for every t S [0, +cxd) p" converges to pt in CPp(M), weakly converges to ptVt 
in M(M); moreover, p^ converges to {pt,ptVt) = '5^t{tA '\^p(I^) for every t € [0, +oo) \1{p). 

(f) For every < s <t there exists a ps-essentially unique monotone map x* g Lp^(R) such that 

(2.18) = (x*)#p3, lim^^— mL2(R), = x, 

(2.19) vt{y)= f vAx)dpl^\x) = it~s)-'(y- f xlix)dp;-^'ix)) for pt-a.e. y eR, 

where Py~*^ is the disintegration of ps with respect to x* . 

Let us recall that the disintegration Py^*" of ps with respect to the Borel (monotone) map 
X* is a Borcl family of parametrized measures uniquely determined for pt-a.c. y £ R, such that 
= ImPr^'^Ptiy) ^ith Pr*((><^)"'(y)) = 1 (see e.g. [1, Thm. 5.3.1]). 

Notice that for a fixed t the map J^t '■ ^p{^) — > '\^p(R) may fail to be continuous with respect 
to the distance Dp, at least in the momentum component pv. 

The gradient flow of the (opposite) squared Wasserstein distance. (2.18) and (2.19) show 
an interesting connection between the semigroup S't in "^2 (R) and the gradient flow "^^f in J'2 (R) 
of the (opposite) squared distance functional 

(2.20) <P^ip)-=-lwi{p,a) Vp,<je92{R). 

Let us recall [1] that for every choice of a reference measure cr G CP2(M) it is possible to define 
a unique continuous and 1-expansive semigroup : CP2(]R) — > CP2(R), t > 0, whose Lipschitz 
trajectories pr ■— ^!^{p) can be uniquely characterized by the Evolution Variational Inequality 

(2.21) ~W^{p,,r^)-]^Wi{p,,r^)<r{ri)-r{Pr) Vyy G T2(K). 
The next result shows that ,5^t and ^/"o basically coincide, up to the rescaling 

(2.22) T = \0gt, t^e\ Pr=Per. 

Theorem 2.4 (Gradient flow of the Wasserstein distance and SPS). Let {pt, ptVt) = =^t(po, Po^^o) € 
V2(K) he the semigroup solution of the sticky-particle system. The Lipschitz curve {pt)t>o in 3'2(R) 
for a.e. t > it solves the Evolution Variational Inequality 

(2.23) i±wi{pu7j)-^Wi{pt,ij)<r"{v)~r°{pt) a.e. zn(0,+oo), Vr; G J'2(R). 
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Equivalently, the reparametrized solutions Pr ~ p^-r satisfy (2.21) with a Pq and we thus get the 
representation formula 

(2.24) pr='^{!lsPs or, equivalently, pt = '^^fj'^^^^^^p^ Vr = log* > 5 = loge. 

Conversely, if t pt is a Lipschitz curve in 'J'2(M) satisfying (2.23) and the initial velocity 
condition 

(2.25) limt-^ \x + tvo{x)-y\^dpt{x,y)=0 Pt = To{po, pt), 

Jr 

then there exists a unique Borel velocity vector field Vt 6 ip^(M) such that {pt,ptVt) = '5^t{po, PoVq) . 
Vt is the Wasserstein velocity field of pt [1, Thm. 8.4.5]. 

Notice that (2.25) corresponds to (2.18) for s = in the case (which a posteriori is always 
verified) = (i x x^)#po- 

We can use (2.24) to exhibit the sohition pt of SPS by a simple limit procedure: 

Theorem 2.5. Let {pt,ptVt) — '^t{po, PqVq) S V2(IR.) be the solution of SPS and let p^ := (i + 
£Vo)^po, £ > 0. Then 

(2.26) pt=hm^#i^°,,/^)(p,) zn5'2(M). 

Moreover, if for some Eq > the map i + SqVo is p^- essentially nondecreasing then 

(2.27) p, - p„ Pt = Kl(t/e}iPe) Ve e (0,£o], < > e. 

The evolution in Lagrangian coordinates. We conclude this section with an even more ex- 
plicit formula for the evolution of the monotone rearrangement function X{t) = Xp^. We denote 
by Ix the indicator (convex, lower semicontinuous) function of % in L^{0, 1) 

(2.28) Ix{X) = liX eX, ^.^j^ subdiffercntial dix ■ i^(0, 1) ^ 2^'^°'^\ 

I +0O otherwise, 

We also introduce the closed subspacc Jix C L^{0, 1), X G 3C, whose functions Y G L^(0, 1) are 
essentially constant in each open interval (a, b) C (0, 1) where X is constant: it is not difficult to 
check that for every X eX and Y G £^(0, 1) 

(2.29) Ye'Kx iff y = 2/ oX for some Borel map y G ip(M), p = X#A. 

Theorem 2.6 (Lagrangian evolution). A curve {pt,ptVt) G V2(K), t>Q, is the semigroup solution 
S^t{po, PqVq) of SPS as in Theorem 2.2 if and only if its monotone rearrangement X{t) ~ Xp^ G 
X C L^(0, 1) satisfies one of the following three (equivalent) characterizations in terms of the 
couple Xq := Xp^ and Vq := ?;o(-'^o) G '^Xo- 

I. X is the unique strong (i.e. absolutely continuous) solution of the Cauchy problem for the 
subdiffercntial inclusion 

(L.I) G -dIjciX) + Vo, X{0) = Xo. 

II. X admits the representation formula 

(L.II) X{t)^PociXo + tVo) 

where Pjc is the -projection on the convex cone X C L^{0, 1). 
III. X is the unique strong solution of the rescaled gradient flow 

(L.IU) t—X(t) e -dIx(X(t))+X(t]~Xo, suchthat \imt-\X(t) - Xq) = Vq in L'^( 0,1). 
at tio 

In each of these cases the curve t i— > X(t) is Lipschitz continuous in L^{0, 1) and right- differentiable 
at each time t; the velocity field Vt can be recovered by the formula 

(L.a) Vit)^^^X{t)^vtoX{t)^P:K^^^^iVo)e^xit) Vt > 0, 
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where P'Kx denotes the orthogonal projection on the closed subspace 3ix C L^{0, 1). The closed 
subspaces IKx(t) are nonincreasing 

(L.b) 'Kx(t)(^'Kxis) tfO<s<t, 

and V satisfy the semigroup identities 

(L.c) X{t)^Poc{X{s) + {t~s)V{s)), V{t) = P:Kx^,^{V{s)) VO < s < i. 

This result shows that the natural evolution space for the Lagrangian sticky particles flow is 

(2.30) Xp{0,l) -.^ \^{X,V) e LP{0,1) X LP{0,1) : X eX, V = voX e^Hx^ p>2, 
endowed with the product distance in LP{Q, 1) x LP{0, 1). The bijective map 

(2.31) (p, pv) e Vp(R) ^ (X, V) e Xp(0, 1), X = Xp, V^vo Xp, 
is in fact an isometry with respect to Dp of (2.10). 

Corollary 2.7 (Lagrangian semigroup). For every p > 2 the time dependent transformations St '■ 

Xp(0, 1) Xp(0, 1), t > 0, which map a couple (Xq^ Vq) G Xp(0, 1) into the couple {X{t), V{t)) — 

I 

St{Xo, Vq) € Xp(0, 1) where X is the solution of (one of the equivalent) (L.I, II, III) and V = X 
as in (L.a), define a right- continuous semigroup in Xp(0, 1), satisfying 

iXit),V{t))^St{Xo,Vo) ^ ipt,PtVt) = ytiPo,Povo) 
^^■^^^ where pt = {X{t))^X, V{t) = Vt o X{t). 

Remark 2.8 (Rescaling). Up to the rescaling r = logt, X{t) = X{e'^), (L.III) is equivalent to 

(2.33) -^X(t) e -dl3c{X{t)) + X{t) - Xo. 

dr 

We shall show (see Theorem 3.1) the Pjc is a contraction in every LP{0, 1), so that (L.II) provides 
a simple and sharp way to estimate X{t) in terms of the initial data corresponding to (2.14b). 
Applying a general result of [22, 23], one can obtain (2.14c) from the representation (L.I). 

Let us finally remark that the Wasserstein gradient flow of Theorem 2.4 is equivalent to (L.III)- 

(2.33) : it is sufficient to introduce the functional 

(2.34) $-(X) :=-i||X-X,|U2(oa)+/3c(^), XeL2(0,l), 
which is related to (p"' by 

(2.35) r{p)^'^''{Xp) Vpea'2(M), 

and is a smooth quadratic perturbation of the convex and lower semicontinuous indicator functional 
Ijc', since 

(2.36) a$"(X) = dIociX) -{X- X^), 

(2.33) is the subdiffercntial formulation in i^(0, 1) of the gradient flow of whose metric 
characterization [1] yields (2.21) thanks to the isometry p ^ Xp between J'2(]R) and %. 

Remark 2.9 (Minimal Lagrangian description). One can use (as in [25, 2]) the initial measure 
Po G T(M) as a reference for the Lagrangian evolution, thus representing pt as x(i)^po for the 
optimal monotone map x{t) = Xq g L^^(R) according to Theorem 2.3 (f). We can therefore 
introduce the convex set OC{po) of essentially nonincreasing Borel maps in the Hilbert space L^^^ (R) 
and we have the corresponding formulae for the evolution in i^^^ (R) (i : R ^ R denotes the identity 
map) 

(L.F) ^x(t)£-9/3c(po)(x(i))+^;o, x(0) = i, 

(L.IF) x{t) = P3c(po)(i + tvo), Kx) = X, 

(L.IIF) t Ax(t) e (x(t)) + x(t) - i. 
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to be completed with the expression for the velocity ^x(t) = v(t) = Vt o x{t). All these relations 
could be easily deduced by Theorem 2.6, since the correspondence x <-> X = xoXq is an isometry 
between L'^^{R) and the closed subspace JCxq of L^{0, 1). On the other hand, it is easier to deal 
with the convex set % in the space L^{0, 1) with the uniform Lebesgue measure as a reference 
and the description provided by Theorem 2.6 is more general, since it allows to compare solutions 
arising from different initial data. 

3. Main properties of % 

In this section we will study the properties of the convex set 3C of nondecreasing functions in 
L^(0, 1), in particular the L^(0, l)-projection operator Pjc and the subdifferential of the indicator 
function Ijc (2.28). Denoting by (•!•) (resp. || • |j) the usual scalar product (resp. the induced 
norm) in i^(0, 1), since 3C is a convex cone, Pjc can be characterized by 

(3.1) g^Pocif) ^ geT, (/ - g|z - .g) < Vz e 3C 

(3.2) ^ geX, if-g\z)<0 VzeX, if-g\g) = 0. 

The next result provides a useful characterization of P3c(/) hi terms of the convex envelope of the 
primitive of /. Recall that the convex envelope of a given continuous function F : [0, 1] — > M is 
defined as 

(3.3) F**{w) := snp^^a + bw: a,beR, a + bv < F{v) Vwe[0, 1]| it) £ [0,1], 

and it is the greatest bounded, (lower semi-) continuous, and convex function G satisfying G < F 
in [0, 1]; it is therefore right and left diffcrcntiable at every point t £ (0, 1) and its right derivative 
g := F** is nondecreasing and right continuous. 

Theorem 3.1 (Projection on %). Let f E i^(0,l) and let F{w) = Jq .fis)ds be its primitive. 
Then 



Pjc(f)=g = 




where F** is the convex envelope of F defined by (3.3). Moreover, for every convex lower semi- 
continuous function ip : M. {~oo, +oo] and every f,hE L^{Q, 1) we have 

(3.4) [ i^{Pjc{f))dw< f ij{f)dw, f i'{Pjc{f)-Px{h))dw< f iif~h)dw. 

JR JR JR JR 

In particular, 9% is a contraction in every space i^(0, 1), p G [1, +oo\ : 

(3.5) ||P3c(/)-PDv(/i)||Lp(o,i) < II./-.9||l.(o,i) V/,/ieLf(0,l). 
We split the proof in several steps. Here is a preliminary Lemma. 

Lemma 3.2. For every f S L^(0, 1) F** is continuous in [0,1], locally Lipschitz in (0,1), and 
coincides with F at w = and w = 1. If f E L°°{0, 1) then F and F** are Lipschitz continuous 
in the closed interval [0, 1]. 

Proof. Let us first assume / e L°°{0, 1) and let L be the Lipschitz constant of F; then 

F{0) - Lw < F{w), F{1) + L{w -1) < F{w) Vwe[0,l], 
so that F**{0) = F(0), F**{1) = F{1), and 

(3.6) F{Q)-Lw<F**{w), F{1) + L{w - 1) < F** {w) Vwe[0,l]. 

Therefore the right derivative g of F** satisfies -L < g{0) < g{w) < ■^F**{1) < L so that F** 
is a Lipschitz function. 

In the general case when / G L^(0, 1), we can approximate its (absolutely continuous) primitive 
F by an increasing sequence of Lipschitz functions Fn uniformly converging to F, e.g. by setting 

Fn(w) = inf F(v)+n\v — w\. 
ve[o,i] 



10 



LUCA NATILE AND GIUSEPPE SAVARE 



Thus F** is an increasing sequence of Lipschitz functions satisfying F**{w) ~ F^^w) at w = 0, 1, 
and pointwise converging to some lower scmicontinuous convex function G as n f +oo with 

(3.7) G{w) < F**{w) < F{w) Vwe[0,l]. 

On the other hand, for w = 0,1 we have G{w) = Hm„|_|_oo Fn{w) = so that F**{'w) = F{w). 

(3.7) also yields 

G{w) < liminf F**(t>) < limsupi^**(i;) < F{w) Vw G [0, 1] 

so that F** is also continuous at w = 0, 1, where G = F. □ 
Let us now consider the set 

(3.8) A = {u- e [0, 1] : (F - F**){w) > o}; 

since A is open and does not contain and 1, it is the disjoint union of a (at most countable) 
collection of open intervals. 

Lemma 3.3. // (a, &) e is a connected component of A then for every w = (I — 6)a + Oh, 

e [0, 1] 

(3.9) F**{{l-e)a + eh) = {l-6)F{a) + eF{h) Ve'e[0,l], F{a)^F**{a), F{b) = F**{b). 

Proof. Since a,b ^ A one has F{a) — F**{a) and F{b) — F**{b). Let w £ [a, b] a minimizer of the 
continuous function 

w^F{w)~L{w), L{w) :^ F{a) + {w ~ a) ^^^] ~ '^^"^ , 

— a 

so that F{w) > F{'w) + L{w — w) for every w G [a, b]. The continuous function 

^F**{w) ilw(^[a,b], 
max [F**{w), F{w) + L(w - w)) if w e [a, b], 

provides a convex lower bound of F and therefore G{w) < F**{w) for every w G [0,1]. Since 
G{w) = F{w) we deduce that w ^ A and therefore H) coincide with a or 6 and the inequality 
G{w) < F**{w) yields F**((l - 0)a + Ob) > (1 - 0)F{a) + eF{b); the opposite inequality is a 
consequence of the convexity of F**. □ 

The next lemma contains the crucial inequality we need to characterize Pjc- 

Lemma 3.4. Let ip G C"'(M) be a convex function. For every f G i^(0,l) and z e % with 
g {F**)', if if - g)i;'iz - g) G L\0, 1) we have 

(3.11) / {f{w)~g{w))^'{z{w)-g{w))dw<0< f {f{w)~g{w))^'{g{w)-z{w))dw. 
Jo Jo 

Proof. We decompose [0, 1] in the disjoint union of the open intervals (a, 6) G covering A (see 
(3.8)) and of [0, 1] \ A, where F{w) = F**{w), and therefore f{w) = g{w) up to a ^^-negligible 
set (recall that F** is locally Lipschitz). In each (a, 6) G F** is linear, g is constant, and the 
function w ip'{z{w) — g) is bounded and nondecreasing, thus its distributional derivative is a 
nonnegative finite measure 7a. Since F ~ F** in {a, 6}, we have 

(/ - g) ^'{z ~g)dw= f{f- g)i''{z - g) dw + [ (/ - g)4,'{z - g) dw 

JA J[0.1]\A 

The second inequality of (3.11) can be simply obtained by considering the convex function ipir) :~ 
ipi-r). □ 



(3.10) G{w) 
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End of the proof of Theorem 3.1. Concerning the projection in Lp'{0, 1), by a standard approxi- 
mation argument, it is not restrictive to assume / e L°°(0, 1) so that g G L°°(0, 1) too. Choosing 
V'(r) := (3.11) yields (3.1). 

In order to prove (3.4), a standard approximation of ip by the increasing sequence of its Moreau- 
Yosida approximations tpnir) minggR ■0(s) + — shows that it is not restrictive to assume 
ip convex, C^, and at most quadratically growing as |r| — > oo. We can then apply the standard 
convexity inequality tp{s) — ^{r) > ^'{r){s — r) and Lemma 3.4 obtaining 

/ U{f - h) ~ ^p{P:K{f) - Px{h))) dw 

> I V''(P3c(./) - P3c{h)) ((./ - P3c(/)) -{h- Poc{h))) dw ^T' 0. 
The first inequality of (3.4) is a particular case of the second one, with h = Poc{h) =0. □ 

The following result is a simple consequence of Theorem 3.1. Let us first introduce for a given 
/ G i^(0, 1) the open set fi/ C (0, 1) where / is locally constant 

(3.12) ^/ £ (0; 1) • / is essentially constant in a neighborhood of ui}. 
Equivalently 17/ is the complement of the support of the distributional derivative of /. 
Corollary 3.5. Let f G L'^{0, 1) and g = P3c(/)- Then 

(3.13) %cr2g. 

Proof Notice that A C (A has been defined by (3.8)); if w G 17/ \ A then F{w) = F**{w), so 
that any linear part of the graph of F in an open interval containing w should locally coincide 
with F**; it follows that F** = F in a neighborhood of w so that w € fig. □ 

Definition 3.6 (The polar cone and the subdifferential of the indicator function Ix). We denote 
by 3C° the polar cone of %, defined by 

(3.14) /g3C° ^ (/|2)<0 Vzg3C ^ P3c(/) = 0. 

The subdifferential dl-x{g) of the indicator function of % (see (2.34)) at some function g G 3C is 
the subset of L^(0, 1) characterized by 

(3.15) ^edlocig) ^ i£.\z-g)<0 VzeX. 
Remark 3.7. X° and dix are clearly hnked by X° = dlx{0) and 

(3.16) ^edljcig) ^ ie%°, (CI5) = 0. 

3C° provides an equivalent reformulation of (3.1), since 

(3.17) g^Pxif) ^ geK, f-ge%°, {f~9\g) = ^ f-gedlxig)- 

If Vl is an open subset of (0, 1), we denote by J\fn the convex cone 

(3.18) 3Vo := |f G C°([0,1]) : F > in [0,1], F = in [0, 1] \ f^j. 
We can give a useful characterization of %° in term of the cone 3Nf := ^(0,1) • 

Proposition 3.8 (A characterization of the polar cone %°). A function f belongs to the polar 
cone 3C° if and only if its primitive F{w) :— J^^ /(s) ds belongs to 3\f. 

Proof. If F G ?^ then one easily gets for every z G 3C n C^([0, 1]) 

(3.19) {f\z)= ( F'[w)z{w)dw = - ( F{w)z'[w)dw 

Jq Jo 

since F, z' > 0, F(0) = F(l) = 0. 
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Let US now assume that / G 3C° ; for every continuous and noimegative function z > and c S R 
with Z(w) — Jg z{s) ds — c, since Z £ OC we have 

Q>{f\Z)^ [ f{w)Z{w)dw = - [ F{w)z{w)dw + F{l){Z{l)-c) 
Jo Jo 

Since c, z are arbitrary, we conclude that F E "N. □ 

The last result of this section concerns a precise characterization of dljc ■ Let us first define for 
/ G L^{Q,1) the closed subspace Oif C L^{Q, 1) defined as 

(3.20) 3-Cf := {/i G L^{0, 1) : /i is essentially constant in each connected component of fij^}. 
We denote by PjCy the orthogonal L^-projection onJif. It is easy to check that 
PM,(/) = /a.e. in (0,l)\r!g, 

(3.21) 

Po<g{f) = J f{w)dw a.e. in every connected component (a,/3) C fig. 

J a 

Moreover, denoting by F the primitive function of /, 

(3.22) if G then / is orthogonal to 

since f ~ F' vanishes a.e. outside Qg and for every connected component (a, /?) of Qg we have 
jj^f{w)dw^F{P)-Fia)^0. 

Theorem 3.9 (The subdifferential of /k)- Let g G 3C, ^ G i^(0, 1), andE{w) := ^{s)ds. Then 
we have 

(3.23) ^edlocig) ^ EeT^n,- 
In particular, 

^ = a.e. in [0, 1] \ fig, 

(3.24) if ( e diocig) then { 



S^{w)dw = for every connected component {a, (3) of 



9' 



so that ^ is orthogonal to "Kg and we have by (3.17) and (3.13) 

(3.25) g^Pocif) => ,g=PM,(/), MgC^Kf. 

Proof. The left implication in (3.23) is immediate, since S G Jvfn implies S G ?vf and therefore 
^ G 3C° by Proposition 3.8; moreover, ^ is orthogonal to 3-Cg by (3.22) and therefore it is also 
orthogonal to g G JCg, so that ^ G dl-x{g) by (3.16). 

Conversely, if G dljcig), then S G 3^ by (3.16) and Proposition 3.8. Moreover, denoting by 
J = g' the nonnegative Radon measure associated to the distributional derivative of g in (0, 1), 
the next Lemma 3.10 yields 

(3.26) f i{w)g{w) dw ^^=^' - f E{w) d-i{w), 

Jo Jo 

which shows that S(u;) = on the support of 7 and yields S G 3Nfog. □ 

Lemma 3.10. Let g G 3C and ^ G 0C° with (nonnegative) primitive S G If "f = g' is the 
nonnegative Radon measure associate to the distributional derivative of g in (0, 1) then S G ^^(7) 
and 

(3.27) [ g{w)^{w) dw = - [ S{w)d'y{w). 
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Proof. Since 7 is a nonncgatativc Radon measure in (0, 1) but not necessarily finite, we need 
an approximation argument to justify (3.27). Let ipn G C^(0,1) be an increasing sequence of 
nonnegative functions such that Vnnni+oc (pn{w) = 1, \^'n\ ^ 2n, and fn{w) = 1 for l/n < w < 
1 — 1/n. We have 

(3.28a) / g^ipn dw = - Eipn d-y ~ '^gf'n dw 



(3.28b) = - / Elfin d7 - / Egifi'^ dw - Egifi'„ dw 

Jo Jo Jl-l/n 



Applying Hardy inequality, we get 



rl/n 

/ ^9Vn dw 
Jo 



< 



2n\\w ^S||i2(o4/„)||wg||i,2(oa/„) < 2C||^||i,2(oa) ||.g||L2(o,i/„) 



so that the integral vanishes as n | +00. A similar argument holds for the last integral of (3.28b). 
Passing to the limit in (3.28a,b) as n f +00 and using Lebesgue dominated (being g£, G L^{0, 1)) 
or monotone (being S > and ipn increasing) convergence theorem, we conclude. □ 

The last Lemma of this section provides a useful example concerning a class of elements in 

dijcig). 

Lemma 3.11 (An example of minimal selection in dix)- If g, h E OC, then 

(3.29) ^h:^P:KAh)-hedIoc{g). 
Moreover, 

(3.30) \\z-h-^h\\mo,i)<\\z-h-^\\L2^o,i) Vee5/3c(.9), ze^Kg. 
In particular, 

(3.31) ifzeJigthen ||z||l^(o,i) < Ik " eilL^(o.i) G ^/^(g). 

Proof. Since h — PcHg{h) is orthogonal to Jig (thus in particular to g), by (3.16) we have to check 
that G X°, by applying Proposition 3.8. By (3.21), = a.e. in (0,1) \ fig, so that the 
primitive S/i of satisfies 

^hiw) = / 6i(s)ds. 
Jngn(o,w) 

The thesis then follow if we show that for every connected component (a, [3) of fig we have 
S/i(a) = Eh{p) = = min[Q^^] S/j. Since the characteristic function X[o,a) oi (0, a) belongs to 3-Cg, 
we have 

/•a 

^h{a) = / ^h{w)dw ^{h~ PMg(/i)|X(o,a)) - 0. 
Jo 

A similar argument shows that S^(/3) = 0. Moreover, for w G (a, (3) we have 

r (3-21) r 

^h{w) = I (,h{s)ds = (w — a) 4 h{w)dw— / h(w)dw, 

J a J a J a 

which shows that S/j is concave, and therefore nonnegative in (a,/3). 

(3.30) follows immediately by observing that G (!Kg)^ and z — /i — belongs to 3-Cg and 
therefore it is the orthogonal projection oi z — h onto (JCg)^. □ 

4. The Lagrangian formulation of the discrete sticky particle system 

In this section we shall show that the discrete sticky particle system satisfies the three charac- 
terizations of Theorem 2.6 and we prove Theorem 2.2. 

Notation 4.1. Let us recapitulate our basic notation and definitions 

(1) Pi{t) = {mi,Xi{t),Vi{t)), i E I = {!,■■■ ,N}, t > 0, is a solution of the discrete sticky 
particle system; 

(2) the positions of the particles are ordered: xi{t) < X2{t) < . . . < XNit); 



14 



LUCA NATILE AND GIUSEPPE SAVARE 



(3) the sets Ji{t) {j E I : Xj{t) ~ Xi{t)} are nondecreasing with respect to time. They 
correspond to a single particle of mass (t) "^i- 

(4) At each time t we pick up the collection of minimal indexes 

I{t) { min : i = 1, iV} = {ii (<)<... < iN(^t}} C /, 

so that each Ji{t) is of the form e / : ik{t) < j < ik+i{t)} for some k and ('^j (^)) jg/j-j) 
is a partition of /. 

(5) We denote by < ti < <2 < ■ • ■ t/i < • ■ • < tn-i the (finite) sequence of times at which 
the cardinality of some Ji{t) has an increasing jump; setting to = and tn ~ +oo, 
{[tfnth+i)}j^^Q is the associate partition of the positive real line with step sizes 6h '■= 
th — th-i- 

(6) The functions Xi are continuous and piecewise linear on each interval \th,th+i)) with 
piecewise constant, right continuous derivatives Vi{t) satisfying (1.2). Each set Ji{t) and 
I{t) is also constant in each interval [i^,t/i+i). 

Let pt = X]iG7 'm'i^xi(t) be the measure induced by the discrete sticky-particle system. In order 
to write explicitly the function X{t) := Xp^ we consider the subdivision of [0, 1] given by 

i 

(4.1) = < U'l < . . . < Wn = 1, Wi = + mi = rrij, i E I. 
We also set 

(4.2) w:, := w,{t) = [j Wj, z e /, 
and we notice that 

AT ^ AT 

(4.3) X(t) =^x,(<)1m/., -X{t) = V{t)=J2v^{t)lw.. 

i=l i=l 

The main result of this section is 

Theorem 4.2 (Lagrangian formulation of the discrete SPS). The couple {X,V) defined by (4.3) 
satisfies the equations (L.I, II, III) and the properties (L.a,b,c) of Theorem 2.6. In particular, it 
defines a semigroup St in the discrete suhspace 

N 

(4.4) X := \{X,V) e Xp(0, 1) : X = ^Xilw, for a finite interval partition (Wi),^i o/[0, 1)| 

i=l 

We split the proof in various steps. 

The collection (W^i(0)ig/(t) ^ partition of [0, 1). In 1/^(0, 1) we introduce the decreasing family 
of finite dimensional spaces J{(i) whose elements are piecewise constant on each interval Wi{t), 
i G I{t). Notice that, by the very definitions of ^x{t) and 'Kx{t) (3.12) and (3.20) 

(4.5) nx(t)^{OA)\{w^:i(.I{t)}, 'K{t)^%x{t)- 
Besides (4.3), the crucial features describing the evolution of X{t) are 

(4.6) X{t)EXr\'K{t), V{t)E'K{t), %{t)=%h,V{t)^Vh\it(^[th,th+i), 

and the update rule for the velocity (1.2): V{th) is constant in each interval Wi{th) = ^jeJi{th)^j 
and its value is given by 

so that by (3.21) 

(4.7) Vh = P^KAVh-i) ^ P^hAVo) since JCq ^ D 2 D • • • 
which yields (L.a) and (L.b). The next lemma shows (L.III). 
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Lemma 4.3. Let X(t) := Xq + IVq be associated to the free system Pi — {mi,Xi,Vi) given by 
Xiit) = Xi{0) + toi(O), Vi{t) = Vi ^ Vi{Q). Then 

(4.8) X{t) = PM(t)(^(t)) = PM(t)(^o + tFo), 

(4.9) t^^X{t)=tV{t)^X{t)-Xa-E{t) for S(t) -Xq + PM(t)(^o) e a/3c(X(t)). 

Proof. Suppose that t £ [th,th+i); since X{t) £ "Kit) = "Kh C 3<(r) and V{r) = P'H:(r)(Vo) for 
< r < i by (4.7), we have by the hnearity of Pjc(r) 

X{t) = Xo+ f' V{r) dr ^ = ' p^^,, + T l^(r) dr) ^ = ^ PM(t)(^o) + T PM(t)(PM(r)(^o)) dr 



= PM(t)(^o) + / PM(t)(V^))dr = Pj^(t)(Xo + / \/odr) = PM(t)(X(0). 
JO Jo 

From of (4.8) wc have 

I 

t- X{t) = tV(t) * = ^ P^(,)(tyo) ^ = ^ ^Y(t) - P5^(i)(Xo) = X(t) - - S(0, 

where S(t) = Pjf(t)(Xo) — Xq; since Xo G DC and 3i{t) = 3ix(t)i by Lemma 3.11 we conclude that 
S(t) e 9/3c(^(i))- □ 

We conclude now the proof of (L.I) and (L.II); notice that (L.c) follows directly by (L.II) and 
(L.a) via the semigroup property of in ' 



Lemma 4.4. Under the same notation and assumptions as before, we have 

(4.10) U{t) := Vo - V{t) ^Vo- Pjr(t)(K,) e dIjc{X{t)), X{t) = Pk{Xo + tVo). 
Proof. Since {U{t)\X{t)) = (being X{t) e 5£(t)), the first inclusion of (4.10) is equivalent to 

(4.11) U(t) ^Vo- V{t) eX° Vt> 0, 

by (3.16). It is not restrictive to assume that t ~ th and V{t) = Vh for some h £ {1, . . . , H — 1}. 
Since 3<C° is a cone and Vq — Vh can be decomposed into the sum 

h-l 

(4.12) Vo-Vh = J2iVk-Vk+i) 

it is sufficient to prove that Vk — Vk+i G 3C° or, equivalently, that 6k+i{Vk — Vk+i) G 0C° . Since 
Vfc+i = Pjffc+i(yfc) wc obtain 

Sk+iiVk - Vk+i) = 5k+iVk - PM.+i(4+iVi.) = {Xk+i - Xk) - P'K.+i(Xfe+i - Xk) 

= Xk+i - P:K,^,iXk+i) + PM.+,(Xfe) - Xk = P^.^AXk) - Xk G DC" 

by Lemma 3.11. 

The second identity of (4.10) follows now by a similar argument, by checking the conditions 
of (3.17). Since X{t) G X and {X{t) - X{t)\X{t)) = by (4.8), it is sufficient to show that 
X{t) - X(i) G 3C°. On the other hand 

X{t)-X{t)^tVa- [ V{r)dr^ [ {Vo-V{r))dr = [ U{r)dr 
Jo Jo Jo 

and (4.11) shows that U{r) G 3C° for every r > 0. Being 3C° a cone, we conclude. □ 

Proof of Theorem 2.2. Let us now consider two discrete Lagrangian solutions (X^(t), V^{t)) = 
St(XQ, V^^) G X, = 1,2. (3.4), (3.5), and (L.II) immediately yield the estimates 



(4.13) / ^{X\t) - X\t)) dw < / V(Xi - X^ + t{V„' - Vo^)) dw 
Jo Jo 

(4.14) \\X\t) - X\t)h.io,i) < \\X^~X^h.^o,i)+t\\Vo'-V;^h.io,i), 
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which are equivalent to (2.14a) and (2.14b). (L.I) yields ([22, Theorem 3], [23, Theorem 1.2]) 
(4.15) /* \\V' - V^^dr <C{l + t)[Y. Il^oll + 11^0 ll) (ll^o - ^oll + W,' - 

1=1,2 

which is equivalent to (2.14c). □ 

5. Stability and uniqueness of Lagrangian solutions 

Our first result concerns the stability of Lagrangian solutions to (L.I, II, III) of Theorem 2.6 
(in particular it applies to those obtained by the discrete SPS in X). 

I 

Lemma 5.1. Let X",^" := ^ curves satisfying all the equations (L.I, II, III) and the prop- 
erties (L.a,b,c) stated in Theorem 2.6 with respect to initial data , — Vq{Xq) converging to 
X„,Vo^v„{Xo) m LP{0,l),p>2. 

(a) X"{t) converges to X(t) in i^(0, 1), uniformly in each compact interval; X is Lipschitz con- 
tinuous with values in L^(0, 1). 

(b) The Lipschitz curve X is right- differentiable at each point t, with right- continuous derivative 
V{t), and it satisfies (L.I, II, III) and (L.a,b,c) of Theorem 2.6. 

(c) V"" strongly converges to V in L^(0,r;L^(0, 1)) for every T > 0. 

(d) The curve X is differentiable in LP{0, 1) and V is continuous at each point of (0, +oo) \ T, 
where T is the jump set of the nonincreasing map t i— > ||T^(i)||L2(o,i)- 

(e) If V is any weak accumulation point ofV"'(t) in LP{0,1), then Pjfx(t)(^) ~ ^(^)- 

(f) ^"(0 V{t) in LP{0, 1) for every t E [0, +oo) \ 7. 

Proof, (a) is an immediate consequence of (L.II) and (3.5), which also show that X" is uniformly 
Lipschitz continuous with values in LP{0, 1) and Lipschitz constant bounded by || Vq" ||/,p(o.i) . The 
convergence is therefore uniform in each compact interval and the limit function X satisfies the 
same Lipschitz bound with constant ||Vb||LP(o,i). 

(b,c) Standard stability results for gradient flows in Hilbert spaces [10] show that X solves (L.I) 
and (L.III); in particular X is right differentiable in L^(0,1) at each t > 0, with L^(0, l)-right 
derivative V{t) which is rig ht-continuous. (4.15) shows that V is the hmit of K in L^{0, T; L'^{0, 1)) 
for every T > (this proves point (c)): in particular, up to the extraction of a suitable subsequence 
Uk, we can find an .if -^-negligible set N C (0, +00) such that Vn;,{t) V{t) in L^(0, 1) for every 
t € [0, +00) \ iV as fc t -|-c». Passing to the limit in (L.c) and in (L.I) we obtain that 

(5.1) X{t) = Poc{X{s) + {t^s)V{s)), ^X{t) = Vit)E~dIxiXit)) + Vis) 

at 

for every s G [0, +00) \ N and t > s. Since V is right continuous, (5.1) eventually holds for every 
< s < i. 

The projection formula of (5.1) shows that 

(5.2) \\X{t + h)-X{t)\\L.io,i)<h\\VmL.io,i)<h\\V{s)\\L.io,i) VO < s < h>0, 
and, more generally, 

(5.3) [ i^{h-\X{t + h) - X(t)))dw < [ i^iVs)dw<[ ^/'(Vb)dw yO<s<t,h>0 
Jq Jo Jo 

for every convex nonnegative function t/' : M ^ M. (5.2) and the right-differentiability of X in 
L^{0, 1) yields that V{t) is also the right derivative of X in £^(0, 1), its LP norm is not increasing, 
and by (5.3) the family Vg is uniformly p-integrable (by Dunford-Pettis criterion, it is sufficient to 
choose a convex function -0 with ip{r)/\r\P +00 as \r\ +00 and ip oVq G L^{0, 1), see e.g. [21, 
Lemma 3.7]) 

From (L.III) we deduce that tV{t) = X{t) — Xq — S(i) where 5(i) is characterized by 

(5.4) E{t)EdIociX{t)), \\X{t)-Xo^E{t)\\<\\X{t)-Xo-a V ^ & dIjc{X{t)). 
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Applying Lemma 3.11 with g :~ X{t) and h := Xq, wc obtain :^(t) = P^^^^^{Xq) — Xq and 
therefore 

(5.5) tV{t)=X{t)-P^,^,^{Xo), V{t)€:Kxit)- 

It foUows by (3.25) that Oixis) 3 Kx{t) if^ < s < i; moreover, by (2.29), there exists a Borel map 
vt e LP^ (M) such that 

(5.6) v{t)^vtox{t), y(t) = Pm,,^,(f(s)) yo<s<t, 

where the last identity follows by the fact that V{t) belongs to JCx(t) and dIoc{X{t)) is orthogonal 
to JCx(i). 

(d) Let 7 be the jump set of the i^-norm of V{t); we show that V is left-continuous at every 
t G (0,+oo) \ T (this also yields the left-diffcrcntiability of X at t). (L.I) provides the minimal 
selection characterization of V 

(5.7) V{t)eVo-dIoc{X{t)), imOllL^(o,i) < IIK)+^||l^(04) yUdIjc{X{t)) Vt > 0. 

Take an arbitrary increasing sequence tn 1 t such that V{tn) — ^ V in i''(0, 1). Since the graph of 
dioc is strongly-weakly closed in L^{0, 1), we have V G Vq — dlK{X{t)). Passing to the limit in 

(5.7) we obtain 

(5.8) ||1/||l2(04) < Jim ||l^(t„)|U2(o,i) = ||V^(t)|lL2(oa) < ||T^o + ClU^co,!) ^ I ^ dI^{X{t)). 

Since dlK{X{t)) is a closed convex set, it follows that V = V{t) and the convergence is strong in 
L^(0, 1) and therefore also in L^(0, 1), since V(tn) is uniformly p-integrable. 

(e) Let nk be an arbitrary subsequence such that V'^''{t) y in i^(0, 1). Passing to the limit 
in the inclusion F"(i) e - dIociX"{t)) we obtain V e Vq - dIx{X{t)). By Theorem 3.9 any 

- (L a) 

element in dIjc{X{t)) is orthogonal to Kx{t) so that PJ{x^^^iV) — Pj£x(t)(^o) " ^(^)- 

(f) Let now t G (0, +oo) \ T and let Uk, V be as in the previous point (e). Up to the extraction 
of a further subsequence (still denoted by Uk), there exists a dense set S C (0, +oo) such that 
V">'{s) V{s) for every s £ 5, so that 

||F|U2(o.i) <limsup||F"-(i)||L^(o.i) <limsup||V^"-(s)|U2(o,i) = \\Vis)\\mo,i) Vse5, s<t. 

k'\-\-oc fct + oo 

Since t is a continuity point for V we obtain by (5.7) 

(5.9) ll^l|L=(oa) < imOllL2(o.i) < l|Vb - C||l2(o,i) S a/x(X(i)), 

which yields^ = V{t), hmsupj,|_|_^ 11^"*° (0 lU^Co.i) ^ II^(OIIl2(o,i): and the strong convergence of 
V"{t) to V{t) in L^(0, 1). The strong convergence in i^(0, 1) follows by the imiform p-integrability 
estimate (5.3). □ 

Corollary 5.2 (Existence of the Lagrangian semigroup). For every initial data {Xq, Vq) G ^2(0, 1) 
there exists a unique Lipschitz curve X in L^(0,1) satisfying the equations (L.I, II, III) and the 

properties (L.a,b,c) stated in Theorem 2.6. Setting V{t) := ■^X(t), the map St : (Xo,Vb) 
{X{t), V{t)) defines a right- continuous semigroup in each space Xp(0, 1), p >2. 

Proof. It is sufficient to approximate {Xq,Vo) G Xp(0, 1) by a sequence {Xq,V^) G X of initial 
data arising from finite discrete distributions of space and velocities in V(R) and to apply the 
previous Lemma. □ 

Corollary 5.3 (Equivalent characterizations). Let {Xq,Vq) G X2(0, 1) be given initial data. If 
X is a solution of one of the equations (L.I), (L.II), (L.III), then it satisfies all the formulations 
(L.I, II, III) and the properties (L.a,b,c) stated in Theorem 2.6. 

Proof. The thesis is obvious in the case of (L.I) and (L.II), whose solution is unique and it should 
coincide with the Lagrangian evolution provided by Corollary 5.2. 

Let us now assume that X is a Lipschitz curve solving (L.III), let X be the Lagrangian solution 
given by the previous Corollary (5.2) with initial data Xq, Vq, and let us set Vjj" := n{X{n^^) — Xo), 
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:= Px{Xq + tVj]"). X"(t) is thus a Lagrangian flow satisfying (L.I, II, III) with respect to 
the initial data Xq, Vq"; in particular 

(5.10) t^^X"{t)e-dIjc{X"it))+X'\t)-Xo, X"(n-i)=X(n-i), 

so that X"(i) = X{t) for t > n^^. On the other hand, the stability Lemma 5.1 yields 

(5.11) \\X'\t)~X{t)\\<t\\V^ ~Vo\\^t\\n{X{n-^)~Xo)-Vo\\ '"^-i"^ as n T +oo, 

so that X ^ X. □ 

6. The continuous sticky particle system in Eulerian coordinates 

In this section we conclude the proofs of the various theorems of Section 2. 

Proof of Proposition 2.1. Starting from (2.9) it is immediate to check that Dp is a metric on 
Vp(]R). Let us check the equivalence characterization (2.13): assuming first that _Dp(/^„,/i) — > 
we obviously have Wp{pm p) ~^ 0; since X„ = Xp^ ^ X ^ Xp and w„(X„) v{X) in £''(0, 1) as 
n I +00, for a continuous and bounded test function C : M R we easily get 

lim / ({x)vn{x) dpn{x) ^ lim / C{Xn {w))vn{Xn{w)) dw 



1 

C{X{w))v{X{w))dw = I C{x)v{x)dp{x), 







(6.1) 

showing that PnVn pv, and 

lim [ \vn{x)\P dpn{x) ^ lim / \vn{Xn{w))\Pdw^ [ \v{X {w))\P dw ^ I \v{x)\P dp{x). 

The converse implication is a particular case of [1, Theorem 5.4.4]: here is a simplified argument. 
If (2.13) holds, then one gets the strong convergence of X„ to X in LP{0, 1); since Vn := Vn o Xn 
is bounded in i^(0, 1), up to the extraction of a suitable subsequence, one has Vn V in LP{0, 1) 
and arguing as in (6.1) 

(6.2) / C{X{w))V{w)dw ^ f C{X{w))v{X{w))dw VC e Cf,(R). 

Jo Jo 

Notice that a function in LP{0, 1) of the form bo X for some Borel map 6 : R R belongs to JCx; 
a simple approximation argument shows that the set o X : C. £ Cb(R)} is dense in JCx so that 

(6.2) yields 

(6.3) voX^P^KxV- 

On the other hand, the last limit property stated in (2.13) yields 

(6.4) ||V^||lp(o.i) < lim ||K||lp(o.i) = IK' o -'^I|lp(o,i) = I|Pmx(^)IUp(o.i) < II^I|lp(o,i), 

nt + oo 

SO that V o X should coincide with V which is also the strong limit of Vn in £''(0, 1). 

Let us finally consider the density of V: if {p, pv) G Vp(R) we can first approximate v in £^(R) 
by a sequence of bounded and continuous functions w„ G Ct,(R). We can then find a sequence 
p^ ~ J2j=i''^j,NSxj N ' N € N, such that p^ p in CPp(R). It is then easy to check that 
VnP^ VnP as t +00 according to (2.13). □ 

Proof of Theorem 2.3. 

(a) The extension of the semigroup ,5^ is not difficult, by using the estimates of Theorem 2.2 and 
the density of V(R) in Vp(R), but not completely trivial since the space Vp(R) is not complete and 
(2.14c)/(4.15) do not provide a pointwise continuous dependence of the velocity from the initial 
data. Therefore, we will use the equivalence stated in Theorem 2.6, which we already proved at 
the level of discrete data in Theorem 4.2, and the Lagrangian stability result of Lemma 5.1. It 
is clear that the only possible extension of S^t to Vp(R) is given by formula (2.32). Since S* is 
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a semigroup in Xp(0, 1) satisfying limtjo St(Xo, Vo) = (Xq,Vo) strongly in LP{0,1)^, S^t satisfies 
(2.f6). 

In order to cheek that S^t is strongly-weakly continuous, we take a sequence /x" = (p", = 
,5^t\[i^\ G with /Xq converging to /i = [p^pv) G Vj,(R) with respect to Dp and we consider the 
associated monotone rearrangement maps {X"-{t),V"{t)) = St(Xo , Vq). By Lemma (5.1) (f), for 
every weakly converging sequence F"* ^ F in L^'(0, 1) and every test function C, G C'°(M) wc have 

Jo Jo 

^'c(X(i))^du;^— '•'^'^^ f\iXit))Vit)dw^=^^ f Cvtdp, 

Jo Jm 

where we used the fact that CiX"{t)) C{X{t)) strongly in LP{0, 1). 

(b) It is immediate to check that (p, pv) = -S^ipoi Po^o) is a distributional solution of (1.1), 
since in Lagrangian coordinates the continuity equation reads 

A C{x) dp,{x) = A £ Q(^x{t)) dw ^ £ C'(^(i))^^(i) 

^"^"^ ^'c'(^(OK(^W)du'= / C'(x)z;t(x)dpt(x), 
and the momentum equation becomes similarly 

±J^Cix)vt{x)dp,{x) ^±j\iX{t))Vit)dw ±j\(Xit))Vodw 

^'c'(X(t))T/(i)Fod«;= /'c'(X(t))t;,(X(i))Fodzi; 
Jo 

(2.29)(L.a) /'\.(^(i))^2(^(^))^^^ f Q' {xV, {x) dpt{x) . 

Jo 

Oleinik entropy condition (1.5) follows easily by (5.5), by observing that P^^^^^^{Xq) is a nonin- 
creasing map, V{t) — vt{X{t)), and pt — {X{t))^X. 

(c) follows from (5.3). 

(d) is equivalent to point (d) of Lemma 5.1; concerning the left continuity of ptVt in the weak 
topology, we fix an arbitrary bounded Lipschitz test function C : M — > M and we observe that 



\im [ C{x)vs{x)dps{x) =lim f C{X{s))V{s) dw = lim [ C{X{t))V{s) dw 
since X{s) X{t) in L^{0, 1) as s t t. On the other hand, since C, o X{t) g ^x{t) we have 

' C{X{t))V{s)dw= [\{X{t))V{t)dw= [\{X{t))vt{X{t))dw^ [ ax)vt{x)dpt{x). 



(e) has already been discussed in point (a)^ except for the convergence at t G (0, +oo)\T, which 
follows from Lemma 5.1 (f). 

(f) (2.18) follows by the projection representation (5.1) and Corollary 3.5. The limit in (2.18) 
can be obtained in Lagrangian coordinate: 



lim 

hiO 



h-\xl+'' - i) - V, dp, = lim / h-\X{s + h)- X{s)) - Vis) 



dw = 



since 1 1-^ X(t) is right differentiable. (2.19) is an immediate consequence of (5.5), which yields 
(t-s)y(t) =X(t)-Pjr,^.,^,(X(s)) VO<s<t □ 

The proof of Theorem 2.6. follows now by applying Lemma 5.1 and its corollaries 5.2, 5.3. □ 
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Proof of Theorem 2.4- (2.23) follows from a simple calculation starting from (L.III): we intro- 
duce the monotone rearrangement Z of the measure 77 G J'2(R) and we observe that M^|(pt, 77) = 
||X(t) — Z|p (we use the usual notation for {X, V) and we denote by || • || the norm in i^(0, 1)). 
We get for some S(t) e dIjc{X{t)) 

I I 

~ W^ipuv) = ~ \\X{t) -Zf = t{X{t)\X{i) Z) {X{t) -Xo- Eit)\X(t) Z) 

^'f (Xit) Xo\Xit) -Z) = \\\X{t) Zf -\\\Z- Xof + \\\X{t) Xof 

Let us consider now the converse implication: if pt satisfies (2.34) then X{t) — Xp^ satisfies (see 
(2.34)) 

(6.5) ~\\X{t)-Zf-^\\X{t)-Zf<^P"{Z)-<i>f»{X{t)) VZe3C, 

which is the equivalent metric formulation [1] of the differential inclusion (L.III). 
Since = {Xo,X{t))#X, (2.25) yields 

(6.6) \imt-^ [ \Xo + tVo-X{t)\^dw^O, 

*i-0 Jo 

i.e. X{t) also satisfies the initial limit condition of (L.III). Therefore, setting V := j^X = v o X, 
by Corollary 5.3 the couple {X{t),V{t)) coincides with the Lagrangian flow S((Xo,Vo) so that 
{pt,ptvt) = yt{po,PoVo). □ 

Proof of Theorem 2. 5. Let us first notice that when i + EqVq is po-essentially nondecreasing, 
(2.27) follows directly from (2.24), since the collision-free motion pt = (i -I- tvQ)^po for t £ [0,eo) 
is a solution of the sticky particle system. 

Let us now consider the general case, setting p^.t ■= '^\og{t/e) ^P<^)- -'^^^ every e > let us consider 
the convex set of bounded Lipschitz functions 

BL{e) := |u G C°-\m) : sup \u\ < e~\Up{u) < (2£)-i| 

and let £ BL{e) be a minimizcr of 

(6.7) me = min ||i;o - u|l = Hi^o - "ell- 

By standard approximation results, lim^io ™c ~ 0, so that Ug converges to wo- 

By the definition of BL{e) the map i -I- eu^ is monotone, and therefore it is the optimal map 
pushing p to /5g = (i -|-eug)^j0o- The sticky particle solution {pe,tT Pe,tVe,t) -S^tiPQi Po'^e) admits 
the representation (2.24) 

SO that, by the exponential rate of expansion of we get 

(6.8) W2{Pe,t,Pe,t) < e-xp{\0g{t/£))W2{pe,Pe) = -W2{Pe,Pe) < t\\vo ~ U^Wl^ im tm^. 

On the other hand, if {pt,ptVt) = '$^t{po, Poi'o), (2.14b) yields 

(6.9) W2ipe,t,Pt) <t\\va -u^\\l2j^) =tm^, so that W2ipt, Pe,t) < '^m^t, 

and concludes the proof of (2.26). □ 

We conclude this section by showing that the representation-convergence Theorem of Brenier 
& Grenier [9] can be easily deduced by our result, in particular by formula (L.II) of Theorem 
2.6. 
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Theorem 6.1 (Brcnicr-Grcnicr). Let vq G C"(R), let pg , N ^ N, be a sequence of discrete 
probability measures supported in a fixed compact interval [— i?, R] and weakly converging to po in 
J'(M), and let p^ be the solution of the discrete SPS with initial data {p^ ,VqPq). For every t >0 
p^ weakly converge to a probability measure pt, whose distribution function Mt{x) :— pt{{~oo,x]), 
t >0, is the unique entropy solution of 

(6.10) dtM + dMM)) = 0, Af(0) = Mq, 
where the flux function A ; [0, 1] — > M is defined by 

(6.11) A{w):= VQ{r)dr, where Vo-.^v^oX^, Xo := Xp„. 

Jo 

Proof. The convergence part foUows by Theorem 2.3 and we can represent Xt := Xp^ by the 
formula Xt = P3c(Xo + tVo) of Theorem 2.6. Introducing the convex primitive functions Ft(w) := 
/™X((r)dr, Theorem 3.1 yields 

(6.12) Ft = (Fo +tA)** so that {Ft)* = (Fq +M)*. 

On the other hand, since the derivative Xt of Ft is the pseudoinverse of Mt (1-9), a standard 
duality result shows that {Ft)* = Gt where Gt{x) = Mt{y) dy, so that 

(6.13) Gt^ {Fo + tA)* ^ {G*+tA)* 

It was already observed by [9, §4] that (6.13) provides the second Hopf formula [3] for the viscosity 
solution of the Hamilton-Jacobi equation 

(6.14) dtG + A{d:,G) ^ in M x (0, +oo), 

and therefore the derivative Alt — dxGt is the entropy solution of (6.10). □ 
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